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NOTATION 

depth  of  submergence  measured  from  the  centerline  of 
the  body 

gravitational  acceleration 
strength  of  point  source 
velocity  vector 

radial  distance  from  source  or  sink 
time 

forward  speed 

rectangular  cartesian  coordinates 

velocity  potential 
wave  height 

regular  wave  height 

local  disturbance 

steady  wave  height 

unsteady  disturbance 

amplitude  of  the  steady  wave  and  unsteady 
disturbance,  respectively 
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ABSTRACT 

The  initial  surface  wave  pattern  due  to  a  submerged  body 
started  from  rest  with  uniform  speed  in  a  perfect  fluid  is 
studied.  The  mean  surface  elevations  may,  in  general,  be  sep¬ 
arated  into  three  parts:  (i)  A  steady  local  disturbance  which 
travels  both  upstream  and  downstream  and  diminishes  as  the  dis¬ 
tance  from  the  body  increases,  (ii)  a  group  of  steady  regular 
waves  which  travel  downstream  from  the  origin  x=x-x  =  0  to 
x  =  1/2  Ut  with  group  velocity  1/2  U,  and  (iii)  a  time-dependent 
cylindrical  disturbance  which  travels  both  upstream  and  down¬ 
stream  am  nlshes  rapidly  as  time  increases.  For  suffi¬ 

ciently  long  xi.ie  t,  the  general  expression  for  the  unsteady  dis¬ 
turbance  can  be  obtained  by  the  method  of  stationary  phase-  It 
is  found  that  for  a  given  finite  value  of  t,  the  effect  of  the 
unsteady  disturbance  becomes  more  serious  as  the  downstream  dis¬ 
tance  from  the  body  Increases,  It  Is  also  shown  that  the  ini¬ 
tial  unsteady  effect  is,  at  lower  Froude  numbers,  more  persist¬ 
ent  and  therefore  takes  a  longer  time  to  subside.  A  comparison 
of  the  predicted  centerline  wave  profile  due  to  a  submerged 
9  feet  long,  7  to  1,  Rankine  ovoid  started  from  rest  with  DTMB 
data  has  been  made;  the  results  are  in  very  good  agreement. 
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INTRODUCTION 

The  problem  of  submerged  bodies  or  ships  in  accelerated 
motion  is  primarily  of  interest  in  connection  with  towing-tank 
experiments.  No  matter  how  quickly  the  final  desired  speed  is 
attained  by  the  model,  the  question  always  arises  as  to  how 
long  it  takes  before  the  effect  of  the  starting  conditions  be¬ 
comes  inappreciable.  This  question  has  been  discussed  by 
Havelock  (l),  (2),  Maruo  (3)  and  Wehausen  (4)  in  dealing  with 
the  wave  resistance  of  submerged  bodies  or  of  thin  snips.  The 
effect  of  the  initial  acceleration  upon  the  wave  pattern  of  sub¬ 
merged  bodies  has  received,  however,  very  little  attention.  It 
is  the  purpose  of  this  report  to  investigate  the  surface  wave 
pattern  of  simple  submerged  bodies  started  from  rest.  The  re¬ 
sult  is  of  practical  importance  in  connection  with  the  experimen¬ 
tal  determination  In  towing  tanks  of  the  surface  waves  made  by 
submarine  models.  The  theoretical  results  for  the  particular 
case  of  a  Ranklne  body  are  compared  with  experiments  conducted 
at  the  David  Taylor  Model  Basin. 

GENERAL  FORMULATION 

The  governing  equation  for  submerged  bodies  moving  through 
an  incompressible,  invlscld  fluid  in  irrotatlonal  motion  may 
be  shown  to  be 


v»<p  =  +  = 

5  x3  dy3  d  z3 


0 


[13 
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where 


x, y,z  are  the  cartesian  coordinates  and  the  origin 
is  taken  on  the  mean  free  surface  with  0-z 
vertically  upwards, 

cp  Is  the  perturbation  velocity  potential  and 

is  assumed  to  be  small. 

The  perturbed  fluid  velocity  q(u,v,w)  is  given  by 


q  =  -  Vcp  =.  -  grad  cp 


[2] 


The  linearized  boundary  conditions  at  the  free  surface  may  be 
expressed  as 


dcp 

at 


g£  =  const.  (dynamical  condition)' 


32.  +  51  = 
az  at 


(kinematic  condition) 


J 


>  at 


0  [33 


where  C  Is  the  free  surface  elevation, 

g  the  gravitational  acceleration,  and 
t  the  time  variable. 

At  the  rigid  submerged  body  surface  it  is  evident  that  the 
condition 


a(cp  +  Ux) 
an 


0 


[4] 
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a 

must  held,  in  which  —  represents  differentiation  along  the  normal 
to  the  boundary  surface  and  where  U  is  the  free  stream  velocity. 

Boundary  condition  £4}  can  generally  be  satisfied  by  dis¬ 
tributing  appropriate  singularities  inside  or  on  the  surface  of 
the  body.  In  order  to  satisfy  boundary  conditions  [3]  and  [4] 
simultaneously,  a  complementary  function,  in  addition  to  the  pre¬ 
scribed  singularities,  is  needed.  It  can  be  shown.  Reference  3> 
that  the  total  velocity  potential  must  have  the  form 


cp  =  cpQ  +  cpj. 


„TT 


H  ( z  —  iu) ) 

=  cp^-  i/  dt  /  d©  /  G(h,0  ,  t  )  e  sin{Vgx(  t-T  )}  Vi*  dK  £  53 


o  TT 


where 


is  the  velocity  potential  due  to  the  pre¬ 
scribed  singularities  and  depends  mainly 
on  the  body  shape 

is  the  image  of  cp 


G(k,9,t)  = 


*0  dn 


k  ( z+lii) 
e  dS 


S' 


uj 


is  an  arbitrary  surface  enclosing  the  singu¬ 
larities 


=  x  cos  0  +  y  sin  9  . 
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Consider  now  the  special  case  of  a  Rankine  body  (ovoid) 
submerged  in  water  at  a  given  depth, as  shown  in  Figure  1.,  sud¬ 
denly  started  from  rest  with  a  given  velocity  U  and  maintaining 
that  speed.  Fcr  a  first  order  approximation, the  fluid  motion 
may  be  taken  to  be  that  due  to  a  point  source  and  a  point  sink 
of  strength  | M |  distributed  at  the  points  (0,0, -f)  and  (C,t,-f) 
respectively.  The  velocity  potential  cpQ  due  to  an  impulsive 
point  source  or  sink  at  (xn*y{  >zn)  l0  given  by 

M 

%  -  r 

n 


where 

M  (t)  »  0  t  s  0 

n v  ' 

M  (t)  =  M  =  constant  t  >  0 
n  '  n 

■W(w, /  +  (y-yn>’  +  (Z-Zn>’ 

The  function  G  of  Equation  [53  in  this  case  la  simply 

h[z  +  i(x  cos  0  +  yn  sin  6)3 

G  =  M  e 


C73 


Substituting  Equation  C 7 3  into  Equation  [ 53  gives 


The  total  velocity  potential  for  a  Rankine  ovoid  started  impul¬ 
sively  from  rest  with  uniform  velocity  can  then  be  shown  to  be 

Cp  =  Cp  +  Cp 

source  ^sink 


[93 


with 

iV  =  x  cos  8  +  y  sin  8  &2'  =  (x+t)cos  8  +  y  sin  8 

l  =  distance  between  the  point  source  and  sink 
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SURFACE  WAVE  PROFILES 


It  is  seen  from  Equation  [3]  that  the  mean  free  surface  eleva¬ 
tions  can  be  expressed  as 


f  ■  1 
D  ^3  J 

g  ot 


z  =  0 


H  (  Z -lUU  )  __ 

xG(x,0,T)e  cos<|-i/gx(t-T  )>  dx 


at  z=0 


[10] 

Consider  first  the  case  of  an  Impulsive  point  source  started  from 
rest  with  velocity  U  which  is  maintained  constant;  the  mean  sur¬ 
face  wave  height  can  then  be  expressed  as 


C  = 


-xf  lxi 

x  e  e  cosVgx(t-T) 


with 


uO  -  j^x-x^.)  +  U  ( t  — t  )j  cos  0  +  y  sin  0 
Rt  denotes  the  real  part  of 
Integrating  with  respect  to  T,  Equation  [ ll]  becomes 


[11] 
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C 


l 


2k  sec  9 


k  -  sec2  0 
U2 


-Kf+lnC (x-x  )cos9+ysin9] 
e 


-xf+lx[ (x-x  )cos9+ysin9] 
k  '  o 

+  I  e 


l(xUCOS0-/gH)t 

e 

1  (  KUCO80-h/gK  )  t 
e 

-  .  .  ’  + 

XCO80  A/-§-  H 

L.  V  U3 

xcose+y^ 

)d*d9 


=  C  +  C .  c 12] 

s  1 

The  limiting  value  of  the  regular  wave  Cr  as  t  becomes  infinite  may 
be  derived  from  t/-e  principle  values  of  the  integrals  in  Equa¬ 
tion  [  12] .  The  steady  local  disturbance  arises  from  the  first 
Integral  in  Equation  [12].  The  steady  surface  wave  pattern  (Cs=Cr+C^) 
due  to  a  point  source  has  been  discussed  in  detail  In  a  previous  re¬ 
port  by  Yim  (5).  The  primary  interest  here  is  to  evaluate  approxi¬ 
mately  the  effect  the  initial  time -dependent  disturbances  have  on 
the  general  wave  profiles. 

For  sufficiently  large  positive  values  of  t.  the  surface  wave 
heights  due  to  the  initial  disturbances  may  be  evaluated  from  the 
following  double  integral  by  applying  the  method  of  stationary  phase 
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7T/2 


r 


«•  j  d9„ 

-ir/2  o 


k  e 


-Kf 


I  ( Kl'i  -  /gHt)  1(kODi  +  v'gKt) 


HCO80 


K  H  COB0 


rd> 


with 


uux  =  [  (; -x  )  +  Ut]  cos  0  +  y  sin  0  =  RA  cos (0-5)  >  0 


R, 


=y(x  -  xo  +  ut)a  + }? 


5X  =  tan 


-1 


x  -  x  +  Ut 
o 


Evaluating  fir  with  respect  to  h.  Equation  [133  has  the  form 


+tt/2 


1  ?1r  J  1 

7irURt|  1 


exp 


gt3f _ gt3  _  T 

4R!3  cos3  (9  -&!  )  ~  \  ¥rx  cos(6 -5x  )  ~  4~ 


r/2 


COB  0  - 


2Ri  cos  (0  -&i  ) 
Ut 


[14 


The  Initial  wave  heights  £  for  sufficiently  large  time  t,  may 
then  be  snowr.  to  be 
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C 


( 


M 

—  R-t 
U 


1 

i 


\ 


exp 

-gty 

f  -  1 

?ta  ‘ 

L^Ri 

‘  4RL  J 

(VF 

=s£ 

Ri 

-l)V 

1  Ri3 
gta 

^cos  5 

2Ri  | 
"  "ut-) 

4r^ 


-2  /2  - 
U 


Ri 


cos  6X  - 


2Ri 

Ut 


sin 


4Ri 


[15] 


This  formula  however  is  not  valid  In  the  neighborhood  of  Rx  =0 
and  Ri.  =  Ut/2  cos  6X  since  in  these  regions  the  non-ci'rcular  part 
of  the  Integrand  in  Equation  [l4]  varies  rather  rapidly  and  the 
method  of  stationary  phase  therefore  fails  to  apply.  The  cor¬ 
rect  evaluation  of  Ci  around  these  points,  which  is  generally  very 
tedious  and  has  little  bearing  on  the  outcome,  will  not  be  pur¬ 
sued  further . 

The  mean  water  free  surface  elevation  due  to  a  Ranklne  ovoid 
started  Impulsively  from  rest  with  constant  velocity  U  can  be 
obtained  by  superposing  the  values  due  to  a  point  source  and  sink 
with  distance  t  apart.  The  resulting  disturbance  may  be  separa¬ 
ted  into  three  parts  as  Havelock  did  in  (1)  for  the  two-dimensional 
case  (see  Figure  2). 

(l)  A  steady  local  disturbance  (C.)  which  traveled  both  up- 
stream  and  downstream  and  diminishes  as  the  distance  from  the 
body  increases. 
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(ii)  A  group  of  steady  regular  waves  (Cr)  which  travel 

downstream  from  the  origin  x-x  =x  =  0  only  to  x  =  -  1/2  Ut  with 

°  _! 

group  velocity  1/2  U,  and  behave  like  x~2  far  downstream 
( 1  «  xA  <  Ut/2X  )  . 

(iii)  A  time-dependent  cylindrical  disturbance  (C^)  which 
travels  both  upstream  and  downstream  and  diminishes  rapidly  as 
time  increases. 

The  first  two  parts  have  been  studied  in  (5).  The  last  part, 
for  sufficiently  long  time  t,  is  simply 


<i* 


2/2  - 

V  (c°s  Bx-  ^ 


^  £t_a 

?  sin  — 


-gt 


4R- 


4Rl 


sin 


\ 


Ut 


R2  cos  da¬ 


's- 


2Ri 


[16] 


where 


Ri  =  -^/(x+Ut)3  +  y3  ,  R2  =-^J{x+t+ Ut)a  +  y3 

respectively . 

61  -  tan'1fe-)  '  6a  ■  tan_1  (x+rnjf) 

In  order  to  assess  quantitatively  the  initial  unsteady  disturbance 
on  the  overall  wave  motion,  simple  numerical  calculations  are 
made  in  the  following  section. 
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NUMERICAL  RESULTS  AND  DISCUSSION 

For  simplicity  only  the  waves  on  the  centerline  of  the  mean 
water  free  surface  due  to  a  Rankine  ovoid  started  impulsively 
from  rest  are  to  be  discussed  in  detail  here.  From  Equation  [l6], 
taking  y  =  0,  the  unsteady  cylindrical  wave  pattern  for 
sufficiently  long  time  takes  the  form 


where 


4F2  1  + 


=  2/2 


1 

Ut 


AY 

ut  i 


1  + 


UF)11 


2x 

Ut 


M 

V. 


is  defined  as  the  amplitude  of 
the  cylindrical  unsteady  dis¬ 
turbance. 
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in  general,  of  0(1), 


F  =  -  Is  the  Froude  number  based  on  depth. 

Vgf 

It  can  be  seen  that  the  amplitude  of  the  unsteady  cylindrical  wave, 
C  ,  for  a  fixed  value  of  x,  diminishes  rapidly  like  1/t  with  in¬ 
creasing  t.  At  a  given  time  t,  the  value  of  decreases  upstream 
but  increases  downstream  as  the  distance  from  the  body  increases. 

It  is,  therefore,  to  be  expected,  that  the  effect  of  the  unsteady 
cylindrical  disturbance  on  the  general  steady  wave  pattern  would 
be  more  serious  far  downstream.  To  illustrate,  it  is  assumed 
that  the  distance  from  the  body  downstream  is  sufficiently  large 
so  that  the  steady  local  disturbance  is  negligibly  small  and  the 
regular  wave  can  be  evaluated  by  the  method  of  stationary  phase, 
i.e.. 


Cs-C<,  +  Cs 


-  c 


cos 


g  1  x-t-C. 


4-  l 
4 


[18] 
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where 


Is  the  amplitude  of  the  st; ady 
wave  profile  far  downstream 

(1  «  xA  <  Ut/2X  ) 


\  =  -  is  the  centerline  wave  length 

g 

From  Equations  [ 17 ]  and  [l83,  the  amplitude  ratio  of  the  unsteady 
disturbance  and  the  steady  wave  is  found  to  be 


=  Idrl 


[19] 


The  non-dlraensional  value  o,f  Z./Z  ,  at  x/X  =  4,  versus  UtA  for 

X  s 

various  Froude  numbers  F  is  shown  in  Figure  3 .  It  can  be  seen 
that  the  value  of  Z ^  diminishes  very  rapidly  with  increase  in  tr 
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In  Figure  4  the  amplitude  ratio  at  TJtA  =  plotted 

against  xA  for  various  F.  It  is  evident  that  the  effect  of  the 
unsteady  disturbance  on  the  steady  wave  profile  becomes  more 
prominent  (since  the  value  of  gets  bigger)  as  the  wave 

moves  further  downstream.  It  is  also  of  interest  to  note  from 
Figures  3  and  4,  that  the  value  of  Cj/Cg  at  a  given  x  and  t  in¬ 
creases  with  decreasing  F.  This  would  seem  to  suggest  the  fact 
that  the  unsteady  disturbance  takes  a  much  longer  time  to  subside 
at  the  lower  Froude  numbers. 

A  preliminary  experimental  study  of  the  Kelvin  wake  produced 
by  submerged  bodies  has  been  made  at  the  David  Taylor  Model  Basin, 
(6).  In  those  tests,  a  9  feet  long,  7  to  1,  Rankine  ovoid  was 
towed  at  several  submergence  depths,  and  wave  height  measurements 
were  made  at  various  points  in.  the  basin.  In  Figure  5  the  center- 
line  wave  pattern  In  a  test  run,  at  U  =  10  ft/sec,  t  =  13.22  sec., 
and  F  =  1,016,  together  with  the  corresponding  calculated  values 
are  shown.  Since  at  distances  very  far  downstream  (>  1/2  Ut)  the 
wave  profile  is  entirely  contributed  by  the  unsteady  disturbance 
as  discussed  in  the  previous  section  and  is  of  no  practical  Inter¬ 
est,  the  comparison  between  the  theoretical  calculations  and  ex¬ 
perimental  data  is,  therefore,  made  only  from  x  =  0  (near  the 
nose  of'  the  ovoid)  extending  downstream  to  the  neighborhood  of 
1/2  Ut.  It  Is  seen  that  the  theory  and  experiment  are  in  quite 
good  agreement.  The  value  of  F,  for  the  particular  case  studied, 
is  rather  large  and  the  effect  of  unsteady  disturbance  on  the 
steady  surface  wave  profile  is,  1  a  general,  small  as  can  be  seer, 
from  Figure  3  or  4 j  the  measured  and  the  steady  surface  profiles 
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are  quite  similar.  At  lower  Froude  numbers,  say  F  =  1/2, however, 
the  unsteady  effect  is  considerably  larger;  it  is  to  be  expected 
then  that  the  true  steady  surface  wave  pattern  may  have  quite  a 
departure  from  the  measured  values  in  this  latter  case. 
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FIGURE  2-SURFACE  WAVE  PATTERN  OF  A  SUBMERGED  BODY  STARTED  IMPULSIVELY  FROM  REST 
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FIGURE  5 -COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  CENTERLINE 
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13  ABSTRACT 

The  initial  surface  wave  pattern  due  to  a  submerged  body  started 
from  rest  with  uniform  speed  in  a  perfect  fluid  Is  studied.  The  mean 
surface  elevations  may,  in  general,  be  separated  into  three  parts: 

(i)  A  steady  local  disturbance  which  travels  both  upstream  and  down¬ 
stream  and  diminishes  as  the  distance  from  the  body  increases,  (il)  a 
group  of  steady  regular  waves  which  travel  downstream  from  the  origin 
x=0  to  x=l/2  Ut  with  group  velocity  1/2  U,  and  (Hi)  a  time-dependent 
cylindrical  disturbance  which  travels  both  upstream  and  downstream  anc 
diminishes  rapidly  as  time  increases.  For  sufficiently  long  time  t, 
the  general  expression  for  the  unsteady  disturbance  can  be  obtained  bj 
the  method  of  stationary  phase.  It  is  found  that  for  a  given  finite 
value  of  t,  the  effect  of  the  unsteady  disturbance  becomes  more  seri¬ 
ous  as  the  downstream  distance  from  the  body  increases.  It  is  also 
shown  that  the  initial  unsteady  effect  is,  at  lower  Froude  numbers, 
more  persistent  and  therefore  takes  a  longer  time  to  subside.  A  com¬ 
parison  of  the  predicted  centerline  wave  profile  due  to  a  submerged 
9  feet  long,  7  to  1,  Rankine  ovoid  started  from  rest  with  DTMB  data 
has  been  made;  the  results  are  in  very  good  agreement. 
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